Proceedings: Evolution of Large Scale Structure - Garching, August 1998 



THE FORMATION AND FRAGMENTATION OF 
PRIMORDIAL MOLECULAR CLOUDS 



Tom Abel 1 ' 2 , Greg Bryan 2 ' 3 and Michael L. Norman 2 ' 4 
1 Max-Planck-Institut fur Astrophysik, D-85740 Garching 
2 LCA, NCSA, University of Illinois, US-61801 Urbana/ Champaign 
3 Massachusetts Institute of Technology, MA, US-02139 Cambridge 
4 Astronomy Department, University of Illinois, Urbana/ Champaign 



ABSTRACT. Many questions in physical cosmology regarding the thermal 
history of the intergalactic medium, chemical enrichment, reionization, etc. 
are thought to be intimately related to the nature and evolution of pregalac- 
tic structure. In particular the efficiency of primordial star formation and the 
primordial IMF are of special interest. We present results from high resolu- 
tion three— dimensional adaptive mesh refinement simulations that follow the 
collapse of primordial molecular clouds and their subsequent fragmentation 
within a cosmologically representative volume. Comoving scales from 128 kpc 
down to 0.5 pc are followed accurately. Dark matter dynamics, hydrodynam- 
ics and all relevant chemical and radiative processes (cooling) are followed 
self-consistently for a cluster normalized CDM structure formation model. 
Primordial molecular clouds with ~ 10 5 solar masses are assembled by merg- 
ers of multiple objects that have formed hydrogen molecules in the gas phase 
with a fractional abundance of < 10 -4 . As the subclumps merge cooling de- 
creases the temperature to ~ 200 Kelvin in multiple "cold pockets" of the 
merger product. Within these cold pockets, quasi-hydrostatically contracting 
cores with masses ~ 100 Mq and number densities > 10 5 cm -3 are found. 
We find that less than 1% of the primordial gas in such small scale structures 
cools and collapses to sufficiently high densities to be available for primor- 
dial star formation. Furthermore, our results indicate that the formation of 
very massive objects, massive black holes, fragmentation of a large fraction 
of baryons into brown dwars or Jupiter size fragments seems, in contrast to 
various claims in the literature, very unlikely. The expected escape fraction 
of UV photons with hu < 11 eV is very small. 



1 Introduction 

Saslaw and Zipoy (1967) realized 
the importance of gas phase H2 
molecule formation in primordial gas 
for the formation of proto-galactic 
objects. Employing this mechanism 
in Jeans unstable clouds, Peebles 
and Dicke (1968) formulated their 
model for the formation of primor- 
dial globular clusters. Further pio- 
neering studies in this subject were 
carried out by Takeda etal. (1969), 
Matsuda etal. (1969), and Hira- 
sawa et al. (1969) who followed in 
detail the gas kinetics in collaps- 
ing objects and studied the possi- 
ble formation of very massive ob- 
jects (VMO's). Hutchins (1976) then 
looked in greater detail at the ef- 
fects of rotation and asked what min- 
imum Jeans mass can be reached 
in a collapsing primordial gas cloud. 
In the 1980's major contributions to 



this field were made by Rees and 
Kashlinsky (1983), Carr etal. (1984), 
and Couchman and Rees (1986), who 
all studied the possible cosmologi- 
cal consequences of population III 
star formation. Their main conclu- 
sion was that for hierarchical struc- 
ture formation scenarios the first ob- 
jects might reheat and reionize the 
universe and thus raise the Jeans 
mass, influencing subsequent struc- 
ture formation quite dramatically. 
Massive Pop III stars would also pre— 
enrich the intergalactic medium with 
metals. 

It is clear that the uncertainties 
mentioned above arise due to the 
inherently multidimensional, nonlin- 
ear, nonequilibrium physics which 
determine the collapse and possi- 
ble fragmentation of gravitationally 
and thermally unstable primordial 
gas clouds. The computational ex- 
pense of solving the network of chem- 
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ical rate equations forced early stud- 
ies to focus on single cell calcula- 
tions (cf. Hirasawa 1969; Hutchins 
1976; Palla et al. 1983; MacLow and 
Shull 1986; Puy etal. 1996; Tegmark 
etal. 1997) adopting simple collapse 
models. Bodenheimer (1986) was the 
first to address the hydrodynamic as- 
pects of the problem, in spherical 
symmetry. Similarly, using a spher- 
ical Lagrangian hydrodynamics code 
and solving the kinetic rate equa- 
tions simultaneously, Haiman, Thoul 
and Loeb (1996) studied the impor- 
tant question of which mass scales 
are able to cool efficiently enough to 
collapse. Unfortunately, the issue of 
fragmentation cannot reliably be ad- 
dressed in such spherically symmet- 
ric models. 

Multi-dimesional studies of first 
structure formation have only re- 
cently become computationally fea- 
sible (Abel 1995, Anninos & Nor- 
man 1996, Zhang etal. 1997, Gnedin 
& Ostriker 1997, Abel etal. 1998a, 
Abel etal. 1998b). Specifically in the 
context of CDM-type structure for- 
mation models we have studied the 
collapse of high-cr density fluctua- 
tions on small mass scales (Abel 
etal. 1998, herafter AANZ98). 

Here we present first results from 
three-dimensional adaptive mesh 
cosmological hydrodynamics follow- 
ing the collapse and fragmentation of 
the very first objects formed in hier- 
archical, CDM-likc models of struc- 
ture formation. 



2 Simulations 

The three dimensional adaptive 
mesh refinement calculations pre- 
sented here use for the hydrodynamic 
portion an algorithm very similar 
to the one described by Berger and 
Collela (1989). The code utilizes an 
adaptive hierarchy of grid patches 
at various levels of resolution. Each 
rectangular grid patch covers some 
region of space in its parent grid 
needing higher resolution, and may 
itself become the parent grid to an 
even higher resolution child grid. 
The general implimentation of AMR 
places no restriction on the number 
of grids at a given level of refinement, 
or the number of levels of refinement. 
Additionally the dark matter is fol- 
lowed with methods similar to the 



ones presented by Couchamn (1991). 
Furthermore, the algorithm of Anni- 
nos etal. (1997) to solve the accu- 
rate time— dependent chemistry and 
cooling model for primordial gas of 
Abel etal. (1997). Detailed descrip- 
tion of the code are given in Bryan 
& Norman (1997,1998), and Norman 

6 Bryan (1997,1998). 

The simulations are initialized at 
redshift 100 with density perturba- 
tions of a sCDM model with Og = 
0.06, h = 0.5, and erg = 0.7. The 
abundances of the 9 chemical species 
(H, H+ , H" , He, He+ , He++ , 
H2 , H^" , e — ) and the temperature 
are initialized as discussed in Anni- 
nos and Norman (1996). After a col- 
lapsing high-cr peaks has been identi- 
fied in a low resolution simulation the 
simulations is reinitialized with mul- 
tiple refinement levels on the Lan- 
grangian volume of the collapsing 
structure. The mass resolution in 
the inital conditions within this Lan- 
grangian volume are 0.53 Mq in the 
gas and 8.96 Mq for the dark mat- 
ter component. The refinement cri- 
terium ensures the local Jeans length 
to be resolved by at least 4 grid 
zones as well as that no cell contains 
more than 4 times its initial mass of 
0.53 Mq. We limit the refinement to 
12 levels within a 64 3 top grid which 
translates to a maximum dynamical 
range of 64 X 2 12 = 262, 144. As we 
will show below the simulation is not 
resolution but physics limited. 

3 Results 

In Figure |l] we show mass 
weighted, spherically averaged quan- 
tities around the densest cell found in 
the simulation at redshift 19.1 of var- 
ious interesting quantities. In panel 
c) we indicate two distinct regions. 
Region I) ranges from the virial ra- 
dius to r coo i ~ 5 pc, the radius 
at which the infalling material has 
cooled down to T ~ 200 K. For 
most of this region the H2 cooling 
time tfj2 is shorter than the free- 
fall time, t ff = [3tt/(32Gp)] 1 / 2 , as 
is illustrated in panel b) of Figure [3]. 
The H2 number fraction rises from 

7 X 10 -6 to 2 X 10 -4 as the free elec- 
tron fraction drops from 2 X 10 — 4 to 
2 X 10" 5 (panel a). 

At r coo i the free-fall time becomes 
smaller than the cooling time. Also 
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Figure 1. Spherically averaged mass 
weighted profiles around the baryon 
density peak at z=19.1. Panel a) 
shows the molcular hydrogen num- 
ber fraction fu 2 an d the free electron 
number fraction x. Panel b) plots 
the H2 cooling time tn 2 > the time 
it takes a sound wave to travel to 
the center, t cross , and the free-fall 
time t ff = [3tt/(32Gp)] 1 / 2 . Panel c) 
and d) show the temperature, bary- 
onic number density and enclosed 
gas mass as a function of radius. 
The vertical dotted line indicates the 
cooling radius (~ 5pc). The virial 
radius of the 5.6 X 10 6 Mq halo is 
106 pc. The cell size on the finest 
grid corresponds to 0.024 pc. Note 
that the simulation box size corre- 
sponds to 6.4 proper kpc. 



a t r coo i the time it takes a sound 
wave to travel to the center, t cross = 
r/c 3 = 7.6 X 10 6 r pc / v / T^ yrs , be- 
comes shorter than the cooling time. 
The cooling time approaches a con- 
stant value at small radii (high densi- 
ties) due to the transition from non- 
LTE to LTE populations of the H 2 
rotational/vibrational states. 

From the time scales one would 
conclude that below r coo i the gas 
should evolve quasi-hydrostatically 
on the cooling time scale. This is, 
however, not strictly true. It turns 
out that centrigufal forces play an 



important role in the collapse. To 
illustrate this Figure |i| shows the 
mass weighted radial profiles of vari- 
ous energies and force terms. From 
the lower panel of Figure ^ it is 
evident that pressure + centrifugal 
forces dominate gravity in the range 
of 0.3 < r/1 pc < 10 yielding a net 
outward force (deceleration). Indeed 
the gas has zero radial velocity at 
T = 0.3 pc, separating the contract- 
ing core of ~ 100 Mq at smaller 
radii. There are multiple possible ori- 
gins for a "centrifugal barrier" as the 
one observed at r = 0.3 pc. The sim- 
plest explanation being angular mo- 
mentum conservation. However, this 
would imply v± X r = const, and the 
centrifugal force, F c , would be pro- 
portional to r~ 3 . However, F c ^ r~ 1 
holds for more than three orders of 
magnitude, as can be seen from Fig- 
ure ^ However, the cooling insta- 
bility and the continous merging of 
small scale structure offer simple al- 
ternative explanations. 



4 Discussion 

Multiple interesting features of 
the collapsing and fragmenting "pri- 
mordial molecular cloud" are iden- 
tified. Most notably is the fact of 
the existence of a contracting quasi- 
hydrostatic core of ~ 100 Mq. 
Within this core the number densi- 
ties increase from 10 5 to 10 8 cm -3 . 
For densities > 10 6 cm - 3, how- 
ever, three— body formation of H2 
becomes the dominant formation 
mechanism transforming all hydro- 
gen into its molecular form (Palla 
etal. 1983). Our chemical reaction 
network does not include this reac- 
tion and the solution cannot be cor- 
rect at r < 0.1 pc. The most inter- 
esting effect of the tree— body reac- 
tion is that it will increase the cool- 
ing rate by a factor ~ 10 3 leading to 
a further dramatic density enhance- 
ment within the core. This will de- 
crease the dynamical timescales to 
100 years effectively "decoupling" 
the evolution of the fragment from 
the evolution of its "host primordial 
molecular cloud". Silk (1983) has ar- 
gued that (due the enhanced cool- 
ing from the 3-body produced H2 
) fragmentation of this core might 
continue until individual fragments 
are opacity limited (i.e. they become 
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Figure 2. Spherically averaged mass weighted profiles around the baryon 
density peak at z=19.1. The upper panel compares local gravitational, ki- 
netic, and thermal energies. The lower panel compares centrifugal ([(*>} 2 m3 — 
(v r )rms]/ r )> pressure (p~ 1 dP/dr), and gravity (GM(r)/r ) forces. They are 
muliplied with rmjj/kg for illustration purposes. Within 0.3 pc gravitation 
dominates leading to a steeper density profile than in the region from 0. 3 to 
10 pc where the sum of centrifugal and pressure forces dominate. 



opaque to their cooling radiation). 
Verifying this possibility will have to 
await yet higher resolution simula- 
tions. 

Clearly, it is the evolution of the 
~ 100 Mq core that could partici- 
pate in population III star formation. 
If it would form an open star cluster 
with 100% efficiency about 6 X 10 63 
UV photons would be liberated dur- 
ing the average life time of mas- 
sive star (~ 5 X 10 7 years). This is 
about hundred times more than the 
~ 4 X 10 el hydrogen atoms within 
the virial radius. However, the aver- 
age recombination time (nfc rec ) _1 ~ 
5 X 10 5 years within the virial radius 
is a factor 100 less than the average 
lifetime of a massive star. Hence, one 
expects only very small or zero UV 
escape fractions for these objects. 

The H2 column density from 
the core to the virial radius is ~ 
10 20 cm -2 . This is much larger than 
the typical 10 14 cm -2 required for 
self— shielding in stationary photodis- 
sociation regions (cf. Bertoldi and 
Draine 1996). Hence also many sub— 
Lyman limit photons will be ab- 
sorbed locally and cannot reach the 
IGM. The host primordial molecu- 
lar cloud is only transparent for pho- 



tons below ~ 11 eV where no Lyman 
Werner Band absorption can occur. 

These optical depth arguments 
will almost certainly become incor- 
rect once supernova explosion alter 
the hydrodynamics and the chem- 
istry and cooling via collisional ion- 
ization, dissociation and shock heat- 
ing. More detailed understanding of 
the role of local feedback will have to 
await new more detailed simulations. 



5 Conclusions 

We have reported first results from 
an ongoing project that studies the 
physics of fragmentation and pri- 
mordial star formation in the cos- 
mological context. The discussed re- 
sults clearly illustrate the advantages 
and power of structured adaptive 
mesh refinement cosmological hydro- 
dynamic methods to cover a wide 
range of mass, length and timescales 
reliably. All findings of AANZ98 are 
confirmed in this study. Among other 
things, these are that 

a significant number fraction of 
hydrogen molecules is only found 
in structures at the intersection of 
filaments 
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only a few percent of the gas in a 
virializcd halo are cooled to T C 

T - 

The improvement of a factor ~ 
1000 in resolution over AANZ98 has 
given new ininsights in the details of 
the fragmentation process and con- 
straints on the possible nature of the 
first structures: 

Only < 1% of the baryons within 
a virialized object can participate 
in population III star formation. 

The formation of super massive 
black holes or very massive ob- 
jects in small halos seem very un- 
likely. 

If the gas were able to frag- 
ment further through 3— body 
and/or opacity limited fragmen- 
tation only a small fraction of the 
gas will be converted into small 
mass objects. 

The escape fraction of UV pho- 
tons above the Lyman limit 
or in the Lyman Werner band 
should be small due to the high 
cloumn densities of HI (Nfii ~ 
10 23 cm" 2 ) and H 2 (N H2 ~ 
10 20 cm -2 ) in the surrounding of 
star forming regions. 

Cautionary Remark: These lat- 
ter conclusions are drawn from one 
simulation and are to be understood 
as preliminary results from an ongo- 
ing investigation. 

Acknowledgments 

Tom Abel acknowledges sup- 
port from NASA grant NAG5-3923 
and useful discussions with Karsten 
Jedamzik, Martin Rees, and Simon 
White. 



References 

Abel, T. 1995, Thesis, University of 
Regensburg, Germany 

Abel, T., Anninos, P., Zhang, Y., 
Norman, M. 1997a, NewA, 2, 181 

Abel, T., Anninos, P., Norman, M., 
Zhang, Y. 1998a, ApJ, in press 

Abel, T., Stebbins, A., Anninos, P., 
Norman, M.L. 1998b, ApJ, in press 

Anninos, P., Norman, M.L. 1996, 
ApJ, 460, 556 



Anninos, P., Zhang, Y., Abel, T., 
Norman, M.L. 1997, NcwA, 2, 209 
Berger, M.J., Collcla, P. 1989, J. 
Comp. Phys., 82, 64 
Bertoldi, F., Drainc, B. T. 1996, 
ApJ, 458, 222 

Bodenheimcr, P.H. 1986, FinaJ Tech- 
nical Report, California Univ. , Santa 
Cruz. 

Bryan, G.L., Norman, M.L. 1997, 
in Computational Astrophysics, eds. 
D.A. Clarke and M. Fall, ASP Con- 
ference #123 

Bryan, G.L., Norman, M.L. 1998, in 
preperation 

Carr, B.J., Bond, J.R., Arnctt, W.D. 
1984, ApJ, 277, 445 
Couchman, H. 1991, ApJL, 368, L23 
Gnedin, N.Y., Ostriker, J. P. 1997, 
ApJ, 486, 581 

Haiman, Z., Thoul, A. A., Locb, A. 

1996, ApJ, 464, 523 

Hirasawa, T. 1969, Progr. Theoret. 
Phys. , 42, 523 

Hutchins, J.B. 1976, ApJ, 205, 103 
Kashlinsky, A., Rees, M.J. 1983, 
MNRAS, 205, 955 

Mac Low, M.-M. & Shull, J.M. 1986, 
ApJ, 302, 585 

Matsuda, T., Sato, H., Takeda, 
H.1969, Progr. Theoret. Phys., 41, 
840 

Norman, M.L., Bryan, G.L. 1997, 
in Workshop on Structured Adaptive 
Mesh Refinement Grid Methods, ed. 
N. Chrisochoidcs 

Norman, M.L., Bryan, G.L. 1998, in 

Numerical Astrophysics 1998, eds. S. 

Miyama & K. Tomisaka 

Palla, F., Salpeter, E.E., Stahler, 

S.W. 1983, ApJ, 271, 632 

Peebles, P.J.E., Dicke, R.H. 1968, 

ApJ, 154, 891 

Puy, D., Signore, M. 1996, A&A, 
305, 371 

Silk, J. 1983, MNRAS, 205, 705 
Takeda, H., Sato, H., Matsuda, T. 
1969, Progr. Theoret. Phys., 41, 840 
Tcgmark, M., Silk, J., Rees, M.J., 
Blanchard, A., Abel, T., Palla, F. 

1997, ApJ, 474, 1 

Zhang, Y., Norman, M.L., Anninos, 
P., & Abel, T. 1997, in S.S. Holt and 
L.G. Mundy, eds., Star formation, 
near and far, AIP Press, New York, 
p329 



